5,972 research outputs found

    Parallel Computation of Finite Element Navier-Stokes codes using MUMPS Solver

    Get PDF
    The study deals with the parallelization of 2D and 3D finite element based Navier-Stokes codes using direct solvers. Development of sparse direct solvers using multifrontal solvers has significantly reduced the computational time of direct solution methods. Although limited by its stringent memory requirements, multifrontal solvers can be computationally efficient. First the performance of MUltifrontal Massively Parallel Solver (MUMPS) is evaluated for both 2D and 3D codes in terms of memory requirements and CPU times. The scalability of both Newton and modified Newton algorithms is tested

    Domain Decomposition Based High Performance Parallel Computing\ud

    Get PDF
    The study deals with the parallelization of finite element based Navier-Stokes codes using domain decomposition and state-ofart sparse direct solvers. There has been significant improvement in the performance of sparse direct solvers. Parallel sparse direct solvers are not found to exhibit good scalability. Hence, the parallelization of sparse direct solvers is done using domain decomposition techniques. A highly efficient sparse direct solver PARDISO is used in this study. The scalability of both Newton and modified Newton algorithms are tested

    Collective non-Abelian instabilities in a melting Color Glass Condensate

    Full text link
    We present first results for 3+1-D simulations of SU(2) Yang-Mills equations for matter expanding into the vacuum after a heavy ion collision. Violations of boost invariance cause a Weibel instability leading soft modes to grow with proper time τ\tau as exp(Γg2μτ)\exp(\Gamma \sqrt{g^2\mu \tau}), where g2μg^2\mu is a scale arising from the saturation of gluons in the nuclear wavefunction. The scale for the growth rate Γ\Gamma is set by a plasmon mass, defined as ωpl=κ0g2μτ\omega_{\rm pl}= \kappa_0 \sqrt{\frac{g^2\mu}{\tau}}, generated dynamically in the collision. We compare the numerical ratio Γ/κ0\Gamma/\kappa_0 to the corresponding value predicted by the Hard Thermal Loop formalism for anisotropic plasmas.Comment: 4 pages, 4 figures, revtex4; v2: typos corrected, discussion on growth rate in expanding system added, accepted for publication in PR
    corecore